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Abstract. Matern's hard-core processes are valuable point process models in spa- 
tial statistics. In order to extend their field of application, Matern's original models 
are generalized here, both as point processes and particle processes. The thinning rule 
uses a distance-dependent probability function, which controls deletion of points close 
together. For this general setting, explicit formulas for first- and second-order charac- 
teristics can be given. Two examples from materials science illustrate the application 
of the models. 
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1 Introduction 

The present paper aims to generalize Matern's well-known first and second hard-core 



point processes. These point process models, introduced by B. Matern Matern, 1960 



1986 , are typical examples for models derived from Poisson point processes, the latter 
being an important basis for constructing more complicated point processes, random 
sets and fibre processes at all Stoyan et al. 1995 . They have been successfully ap- 



plied to real data, for instance, they have been used in ecological [e. g. Picard et al. 



analysis Stoyan 1988 



2005 




Stoyan 


1987 Warren 


1972 and CSMA network modeling [e. g. 


Blaszczyszyn 


2009 


Busson and Chelius 


2009 


Haenggi 


2011 as well as 



Baccclli and 



Matern's first and second hard-core point processes are derived by applying a spe- 
cific thinning rule to a homogeneous Poisson point process in R rf . As such they are 



important examples for dependent thinning Illian et al. 2008 Stoyan 1988 where 
the thinning depends on the underlying process, somehow. For instance, the Matern 
I hard-core point process is obtained by deleting every point in the process with its 
nearest neighbor closer than a given hard-core distance [cf. Matern 1960 19861 . In 



general, a thinning operation or rule Illian et al. 2008] determines which points in 



the basic process are deleted. For example, such thinnings drive the evolution of plant 
communities due to competition-induced mortality |Batista and Maguire 1998| . 

In contrast to depend thinning, the well-known p- and p{x)-thinning approaches 
described in Illian et al. 2008 Section 6.2.1] or Daley and Vere- Jones 2008 Section 



11.3] use independent thinning. That means that the thinning operation is independent 
of the configuration of the underlying point process and, at position x, a point will be 
deleted with some deterministic probability 1 — p or 1 — p(x), respectively. 

In the following, the way of thinning a Poisson point process in order to obtain a 
Matern hard-core point process is generalized in two directions. The first idea is to 
combine both independent and dependent thinning which can simply be interpreted 
as a subsequent independent thinning of a dependently thinned point process. The 
second idea is that in the dependent thinning a distance-dependent probability function 
/ controls deletion of points which are close together. This means that, depending on 
the distance to its neighbors, a point will be deleted only with some probability and 
not surely as is the case, e.g., in the Matern I thinning rule. Thus it seems to be 
justified to speak of a probabilistic thinning rule in what follows. 

An application of appropriate probabilistic thinning rules to Poisson point processes 
then leads to a class of point processes which are generalizations of the Matern hard- 



core point processes. As distinguished from Gibbs point processes Stoyan et al. 1995 



Illian et al. 2008 Daley and Vere-Jones 2008 , explicit formulas for first- and second- 



order characteristics can then still be derived. By means of several examples for the 
probability function / it is shown that soft-core, hard-core as well as aggregative point 
processes can be obtained by this approach which reveals its high flexibility. 

This paper is organized as follows. In Section [2] first the generalization of the 
Matern I hard-core point process based on the thinning of a homogeneous Poisson 
process is introduced and first-order and second-order characteristics are given. Based 
on the idea of Stoyan and Stoyan 1985| and Mansson and Rudemo 2002 to extend 
the original Matern II hard-core model by giving each point a random radius, where in 



Mansson and Rudemol 12002] as a special case also a respective extension of the Matern 



I process is covered, in Section [3] the dependent thinning model is further enhanced to 
marked point processes which is useful to model particle systems. Section [4] contains 
then a discussion of a related generalization of Matern-II-type point processes. Finally, 
in Section [5] the applicability of the new models is illustrated by means of two data 
sets. 
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2 Probabilistic thinning model - Matern I case 



2.1 Model description 



Let $ be a homogeneous Poisson point process in R d with intensity A > [see Daley 
and Vere-Jones 1988 on a probability space (f^J 7 , P), po <E]0, 1], and / : [0,oo[ — > 
[0, 1] be a measurable function. Denote by || • |j Euclidean distance in R d . From $ a 
new model $th i s derived by applying the following probabilistic dependent thinning 
rule to <i>. A point x € $ is retained with probability 



p(x,$)=pol[[l-f(\\x-y\\)] 



(1) 



y e<s> 

y=£x 



independently from deleting or retaining other points of <I>. This means that two 
points a distance r > apart delete each other independently with probability f(r). 
Independently from deleting due to pairwise interaction, each (surviving) point is 
(then) deleted with probability 1 — pq. 

Since the homogeneous Poisson process $ is both stationary and isotropic and the 
thinning rule is independent both from location and direction, the thinned point pro- 
cess <I>th is stationary and isotropic as well. 

In the following we will sometimes write MatI[A,po, /] for the distribution of <5> t h- 

Remark 1. Taking po = 1 and f = 1[o,b] f or some R > 0, the point process <&th 
coincides with the Matern I hard-core process with hard-core distance R since then any 
two points a distance r < R apart delete each other with probability 1, i. e. almost 
surely. 

Remark 2. The thinning rule |Ip could be refined by making the retention probability 
dependent also on three-point or further multi-point configurations, or, more general, 
on any functional which depends both on the point a; e $ and the point pattern An 
example for the latter would be the number nn(x,&) of neighbors of x £ $ within a 
certain distance R > which is, e.g., directly used in the definition of Strauss point 
processes Strauss 1975'. In fact, even the choice f = /ol[o.fi] for some fo G]0,1[ 



results in a retention probability 



p(x^)^ Po (l~f a ) n ^\ 



2.2 First- and second-order characteristics 

As is known for the original Matern I hard-core point process, explicit formulas both 



for its intensity and its pair correlation function Daley and Vere-Jones 1988 Illian 



et al. 



20081 M0ller and Waagepetersenl 2003 Stoyan et al.l 1995 can be stated Daley 



Vere-Joiies[ |1988[ |Illian et al.\ |2008i |Matern[ |1M^|1986[ |Stoyan and Stoyan| 



and 



1985 



Although the definition of the thinned point process <&th is more complicated, 

can be 



arguments similar to that given for instance in Daley and Vere-Jones 1988 



used to derive the subsequent expressions in Theorems [3j and |4j for both the intensity 
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and the pair correlation function of <J> t h- The respective proofs are omitted here due to 
the fact that the model MatI[A,j>o> /] a ^ so appears as a special case of another model 
introduced later in Section [31 

Theorem 3 (Intensity). The intensity Ath of the thinned point process <I>th is 
A t h = Ap exp (-\db d f{r)r d ~ l dr\ , 

where bd denotes the volume of the unit ball in WL d . 

In case the integral in Theorem [3j is infinite, the resulting intensity Ath vanishes. 
That is, <!>th contains almost surely no points since the applied thinning is so strong 
that all points of <i> are removed almost surely. Since this case is uninteresting, in 
what follows we thus consider only those functions / which satisfy the integrability 
condition 



oo 

r d ~ 



~ l f{r)dr < oo. (2) 
Let /©/ denote the radial self-convolution of /, i. e. 







[/©/]«= / f(\\x\\)f(\\x-r-v\\)dx (3) 

for r > and v £ R d with ||u|| = 1, which is the d-dimensional convolution of /(|| ■ ||) 
with itself at point r ■ v. 

Theorem 4 (Pair correlation function). The pair correlation function g t h of the 
thinned point process <&th is 

5th = (l-/) 2 exp(A/0/). 

The following two examples illustrate for dimension d = 2 how certain choices of / 
influence the second-order behavior of the resulting thinned point processes. 

Example 5. Let A > be arbitrary and po = 1. Consider MatI[A, 1, f a ,R\ with 

<r <a 
otherwise 



faM r ) 




for R > and a € [0, R] . Then it is plain to see that for any A and R fixed the 
intensity of(fr t h, which is distributed according to MatI[A, 1, f a ,R\, is 

Ath = Aexp(— n\R 2 ), 

i. e., it does not depend on a. However, depending on the parameter a, the pair cor- 
relation function o/^th shows a certain range of second- order behavior, see Figure^ 
For a = the thinning generates a pure soft-core point process, i. e. thinning is the 
stronger the closer point pairs of the initial Poisson process are but each pair distance 
has still non-vanishing probability. In the other direction, the Matern I hard-core model 
is included as the limit for a — > R. A mixture of both hard- and soft-core type behavior 
can be achieved with a €]0, 1[. 
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Figure 1: Left, plot of f a i for a = (dotted), a = 0.75 (dashed) and a = 1 (solid). 
Right, corresponding pair correlation functions for A = 1/2. 
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Figure 2: Left, plot of f a for a — 8 (dotted), a — 2 (dashed) and a = 0.5 (solid). Right, 
corresponding pair correlation functions for A = 2. 



Example 6. Another type of family MatI[A, 1, h a ], where again A > is arbitrary and 
Po = 1, is given by setting 

K(r) = r ° , r>0 

r(i + §) 

/or a € [0, oo[, where T denotes the T -function. Again, for A fixed, all resulting thinned 
point processes have the same intensity 

A th = Aexp(-7rA), 

i. e. independent from a > 0. Here, the thinning based on h a results in aggregation or 
cluster-like processes which is indicated by pair correlation functions with its maximum 
in the origin, see Figure [J[ since points close together are deleted with relatively low 
probability. 

Examples [5] and [6] show that the behavior of the pair correlation function of the 
thinned process <l>th might be adjusted using appropriate functions / such that it is 
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possible to model soft-core, hard-core as well as aggregative point processes [cf. |Illian| 



et al. 2008 . Furthermore, a 'mixture' of these types can be obtained combining the 
corresponding functions /. In summary, this reveals high flexibility of this approach 
with respect to second-order properties. 



3 Probabilistic thinning for marked Poisson processes 



Marked point processes as generalizations of usual point processes are highly relevant 
in practical applications. Many ecological and environmental systems can be described 
by marked point processes [see Gavrikov and Stoyan 1995 . Furthermore, there is a 



large literature on processes of non-overlapping grains in physics and chemistry [see 
Andersson et al.l 2006 Mansson and Rudemo 2002 . Here, Matern hard-core processes 



are equipped with random radii as marks. While in the previous section a dependent 
thinning model generalizing the usual Matern I hard-core point process was introduced, 
the aim of the present section is to carry over this approach to respective marked 
point processes. Again, intensity and pair correlation function of the corresponding 
unmarked point process can be given explicitly. 

Let be a homogeneous independently marked Poisson point process in R rf with 
intensity A and independent and identically distributed (i. i. d.) real- valued marks with 
H as its mark distribution. Furthermore, let po £ ]0, 1], and / : [0, oo[xR 2 — > [0, 1] be 
a fixed measurable function satisfying /(•, m, n) — /(•, n, m) for all m, n G R. 

From 'J a new model ^th is derived by applying the following probabilistic dependent 
thinning rule to 'J. The marked point (x,m) S ^ is retained with probability 



p(x, m, ^) = po J| I 1 ~ f(W x ~ v\\> m > n )] 



(4) 



(i/,n)e* 



independently from deleting or retaining other marked points of ^ . This means that 
two points of distance r > apart with marks m and n delete each other independently 
with probability f(r, m, n), and, again, independently from deleting due to pairwise 
interaction, each surviving point is then additionally deleted with probability 1 — p a . 

Since the thinning rule is again independent both from location and direction, the 
point process \&th of unmarked points of ^th inherits both stationarity and isotropy 
from the homogeneous Poisson process of unmarked points of , f. 

In the following we will sometimes write MatI[A, p,,po, f] for the distribution of \I>th- 
Of course, if / does not depend on the marks the unmarked point process coincides 



with the model MatI[A,poj/] introduced in Section 2.1 and all formulas given there 
appear as particular cases of the subsequent results. 

Example 7. Taking f(r,m,n) — l[o,m+n]( r ); r,m,n > 0, and random marks uni- 
formly from [0,1], we obtain a hard-core point process with random hard-core radii, 
i. e., we can interpret the retained marked point (x,m) as a ball of radius m centered 
in x, see Figure [3| This example was studied by Mansson and Rudemo ^20021 as a 
model for systems of varying- sized, non-overlapping spherical grains. 



G 



Figure 3: Sample of ^th in dimension d — 2 for A = 100, po — 1, /(•, m, n) = 1 [o,m+r. 
and fi the uniform distribution on [0, 1]. 



Theorem 8 (Intensity). The intensity Ath of the thinned point process ^th ?s 
A t h = Ap / exp[-Xdb d / /(r,m, Z)r d_1 dr,u(dZ) ) /x(dm). 



R ^0 



As an example, consider once more the case po = 1 and f(-,m,l) = l[o, m +i] with 
positive- valued mark distribution /i. Then the result in Theorem [8] reads 

/■OO / />00 \ 

Ath = A y exp ^-A& d y (Z + m) rf (i(dl)J (j,(dm) , 



which coincides with the formula in Mansson and Rudemo[ 2002[ Section 2.2, Example 
2.1]. 

Theorem 9 (Pair correlation function). The pair correlation function of the thinned 
point process ^th is 

\ 2 p 2 r r 

9th{r) = —■ to — / / (1 - f(r, m,n)) 2 q m q n q„ un (r) fi(dm) n(dn), r > 0, (5) 
A th JrJb, 

where 



g m , n ( r ) = exp ( A / f(-,m,l)Qf(-,n,l)(r) n(dl) ) , r > (6) 



and 



9m = exp [ -A / / /(||x||,m,Z)dx/i(dZ) 

JK JR d 



(7) 



Note that g t h (r) actually does not depend on p since it appears also as a factor in 
A t h and thus cancels out. 
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4 Probabilistic thinning model - Matern II case 



In the previous sections, a certain kind of dependent thinning has been introduced 
where two competing points of an underlying Poisson process are both removed with 
some probability depending on a deterministic function /. While this is a general- 
ization of the classical Matern I hard-core point process, in the present section we 
aim to generalize the classical Matern II model in a similar fashion. For that purpose, 
weights will be assigned once to all the points and in a competition between two points 
only one of them, namely the point equipped with weight greater than or equal to the 
weight of the other point will be removed with some probability. Again, we can derive 
expressions for both the intensity and the pair correlation function of the resulting 
process of unmarked points. 

Let II be a homogeneous independently marked Poisson point process in R d with 
intensity A and i. i. d. Revalued marks. The first component m of any such random 
bivariate mark (m, v) has distribution \x and the second component v has distribu- 
tion v m which might depend on m. The mark m plays the same role as in Section [3] 
and is sometimes referred to as 'mark' whereas v serves as a weight used in the thin- 
ning procedure and is thus sometimes referred to as 'weight'. Furthermore, let again 
Po €]0, 1], and / : [0,oo[xR 2 — > [0,1] be a fixed measurable mapping satisfying 
/(•, m, n) — /(•, ro, m) for all m, n £ R. 

From II a new model n t h is derived by applying the following probabilistic dependent 
thinning rule. The point (x,m,v) € H is kept as a point of n t h with probability 

p(x, m,v,U) = po Yl [ l - 1 i v > w ) f(\\ x - y\\i m > n )] ( 8 ) 

(y ,n,w)£Tl 
yjix 

independently from deleting or retaining other points of II, where 1{^4} is the indicator 
of event A. This means that if two points with marks m and n are a distance r > 
apart then only the point with weight greater than or equal to the weight of the other 
point is deleted by the other point with probability /(r, m,n). Additionally, each 
surviving point is then again independently p -thinned. 

In the following we will denote by MatII[A, /i, (%) m gK,Po, /] the distribution of n t h- 
Note that the meaning of the weights in the thinning rule ^ is here in accordance 
with the meaning of the respective weights used in the definition of the original Matern 
II hard-core point process in most of the literature Illian et al. 2008[ Stoyan et al. 



1995 Stoyan and Stoyan) 1985] (but not Mansson and Rudemo 2002] ), i. e., they have 



to be understood more (biologically) as times of appearance than importance weights, 
and in a competition the lower weight wins. 

Remark 10. Taking /(-,m,n) = l[o,m+n](")> m,n > 0, and marks according to a 
positive mark distribution /i results in a hard-core process with random hard-core radii. 
This case was also studied by Mansson and Rudemo f200^ as an extension of Matern' s 



second hard-core point process to random configurations of non-overlapping spheres. 
The original Matern II point process with hard-core radius R > can be obtained 
using \x = <5r and v m as the uniform distribution on [0, 1]. 
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Remark 11. The particular choice v rn = Si for all marks m leads back to the model 
^th described in Section^ Here, all points would have the same weight such that two 
competing points delete each other independently with the same probability. 

Denote by F T the cumulative distribution function of a probability measure r on R, 
i.e. F T (t) = t(] - oo,t]), t € R. 

For simplicity assume in what follows that the weight distributions i>i are all (ab- 
solutely) continuous (w.r.t. lebesgue measure) but may still depend on mark I. The 
main effect is that the event l{v > w} has probability zero and hence only one of 
two competing points is deleted with some probability. In particular, this excludes the 
Matern-I-like processes introduced in the previous sections. 

Furthermore, let us recall that the point process H t h of unmarked points of Ilti, 
inherits both stationarity and isotropy from the homogeneous Poisson process of un- 
marked points of II. 

Theorem 12 (Intensity). The intensity Ath of the thinned point process n t h is 



where 



Ath = Ap / / q m (w) v m (dw) fi(dm) 



q m (w) = exp ( -A / F Ul (w) / f(\\x\\,m,l)dxfj,(dl)). (9) 
Jr. J]R d 



In the special case that the weight distribution i>i does not depend on the mark I, 
i.e., v\ = v for some continuous distribution v, Theorem 1 1 2| simplifies to 

A t h = Ap / / ex P ( -XF v (w) / / /(|M|, m, l)dx fi(dl) ) v(dw) fi(dm) 

= Apo / / expl— Au / / /(||x||,m, l)dxfi(dl) I dufi(dm) 
Jr Jo \ J r ^R d / 

f 1 - exp (-A J R / Rd /(||a;||,m, Qda;/j(dZ)) 

~ P °L J R J Rd f(\\xlmJ)dx^dl) Mdm) ' 



by change of variables. In the case pg = 1 and /(-,m,n) = l[o. m +n] formula (10) 



coincides with the result stated in Mansson and Rudemo 2002 Theorem 3.1] 



Remark 13. $From Equation ^ it is easy to see that, due to F Ul (w) < 1 for all 
w, the intensity of the thinned point process n t h is always greater than the intensity 
of a thinned point process according to MatI[A,/x,po> /] from Section [5] with the same 
parameters. 

Theorem 14 (Pair correlation function). The pair correlation function <? t h of the 
thinned point process Hth is 

ffth( r ) = i2° / (I- f(r,m,n))I r (m,n) fj,(drn) fx(dn) , r > 0, 



\ 2 
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where 

I r (m,n)= / / q m (w)q n (t)q m , n (w,t,r)u m (dw)v n (dt) 

with q m (w) from Equation Q and 

q m ,n{w,t,r) = exp ^A J F Vl (min{w, t}) [/(•, to, l)Of(-, n, l)](r) fi(dl) 

Again gth( r ) does not depend on p . In the special case where v\ = v does not 
depend on I, I r {m,n) can be written as 

I r (m,n) = / / c&g* g TO) „(r) min ^*> dsdt = J(a,b,c) + J(b,a,c) 

where 



o Jo 



T , , , 1 1 / exp(a + 6 + c) exp(6) 

J (a, o, c) — 



&(a + 6 + c) a + c \ a + b + c b 

with q TO and q m _ n (r) according to Equation Q and Equation (|6j|, respectively, and 
a = logq m , 6 = logq„ and c = logg m>n (r). 



5 Applications 



5.1 Fontainebleau sandstone 

The first data set is a point pattern describing the pore network of a sample of 
Fontainebleau sandstone. A visualization is given in F igure |4} A detailed description 



how this point pattern was obtained can be found in Sok et al. 2002] . It has been 
further analyzed in the literature, for instance, iTscheschel and Stoyanl 12003] discuss 



second-order characteristics and a certain Euler-Poincare characteristic connected with 
the data. A standard test of the hypothesis that the pattern is of CSR type (complete 
spatial randomness) jlllian et"aT 2008 results in rejection with a p-value of 0.0002. 



The minimum intcrpoint distance in the pattern is 60.7 ^m, and it is just this hard- 
core distance which leads to a rejection of the CSR hypothesis. Consequently, a hard- 
core point process model seems to be more appropriate for this data. Because of the 
low point density Matern-like point processes are promising. 

The plot in Figure [4] shows the estimated pair correlation function g [see Illian et al 



2008] Section 4.3.3] of the data and the pair correlation function g both of a fitted 
Matern I and a fitted Matern II hard-core point process. 

Taking the minimum interpoint distance of 60.7 fjia as an estimate for the hard-core 
distance R (which is even a maximum likelihood estimate), fitting is here easily done 
by estimating the intensity A of the underlying Poisson process as the only remaining 
unknown parameter by the method of moments. That is, due to 



A/ = Acxp(-A6 3 i? 3 ), A// = 



1 



exp(-A& 3 R 3 ) 
b 3 R 3 : 
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Figure 4: Left, 470 nodes of a network adapted to the pores of a 2 x 2 x 3mm sample of 
Fontaincbleau sandstone. Right, the empirical pcf g (dotted), g for Matern 
I (dashed) and for Matern II (solid line). 



for the intensities of the, respectively, Matern I and Matern II hard-core point process, 
an estimate of A can be obtained by solving for A in the equations A = A/ and A = A//, 
respectively, where A = 39.17 mm -3 is the empirical intensity of the data. However, 
Figure [4] shows clear differences between the respective pair correlation functions in- 
dicating that none of the both Matern hard-core point processes is a good model. 



Since pure hard-core point process models of Matern type seemed to be not appro- 
priate we have then fitted a model MatI[A,po> fR.a.b] from Section 



2.1 



with 



fR,a,b(r) 



i exp 



{r-Rf 



0<r<R, 
otherwise, 



with R — 60.7 /im. Here, as parameter estimate the best possible choice of (\,po, a, b) 
was taken, meant in the sense that under the condition Ath = A the contrast 



A(X,p ,a,b) 



M -5th(r)] 2 dr 



(11) 



is minimized with respect to (\,po, a,b), where [r m i n ,r max ] is a suitable domain. This 
is a variant of the well-known minimum contrast method for parameter estimation 
DiggleJ |2003[ |Heinrictij 1992] where here the difference to be minimized depends on 
the pair correlation function as that summary statistics which is available at least 
via numerical integration for the models under consideration. The minimum contrast 
method using the pair correlation function has been also successfully applied by |Stoyan| 
and Stoyan 1996 and M0ller and Waagepetersen 2003 p. 183]. This resulted in 



estimates A = 1919 mm -3 , po = 0.92, a = 6.3 and b = 3917 /mi 2 . Figure [S] shows the 
estimated function jn.,a,b as well as both the empirical pair correlation function g of the 
network data and the pair correlation function g t h of the estimated model. The visual 



11 




50 100 150 200 250 300 100 200 300 400 



Figure 5: Left, plot of fn } a,b, R = 60.7/zm, a = 6.3 and b — 3917 fim 2 . Right, pair 
correlation functions g (dotted) of the network data and <?th (solid line) of 
the fitted model. 




100 200 300 400 100 200 300 400 

Figure 6: Left, empirical nearest-neighbor distribution function of the network data 
(dotted) and nearest-neighbor distribution function of the fitted model 
(solid). Right, same for empirical empty space function. 



finding from Figure [5] (right) that the fit is good can be suggested by formal tests. 
For instance, a deviation test [see |Illian et al.[|2008[ Section 7.4] for the corresponding 
L- functions Illian et aTj |2008| Section 4.3] with global deviation measure 



\L(r) - L th (r)\ 2 dr 



(p- value 0.18) with k = 99 simulations indicates that now the model was chosen flexible 
enough to give a fit which mimics the second-order behavior of the data sufficiently 
well. Two other deviation tests with k = 99 simulations using the nearest-neighbor 
distance distribution function as well as the empty space function (or 'spherical contact 
distribution function') [Illian et~aL 2008. Section 4.2] (see Figure [6] for plots) instead 
of the L-function then show that the fit is good also in other respects (p-values 0.14 
and 0.12, respectively). 
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Figure 7: One of the used QICPIC images, size 1 m x lm. 



5.2 Patterns of deagglomerated alumina particles 

The second data set are three samples of a mono-layer of deagglomerated alumina 
particles within water which serve as a starting point for the investigation of certain 
agglomeration processes not discussed here. The patterns, one shown in Figure |7j 
were obtained with a QICPIC sensor (Sympatec/Germany), which is a measurement 
device for dynamic picture analysis. For the test setup a liquid dispersing unit was 
used to get such a mono- layer flow of deagglomerated alumina particles through a flat 
cuvette where then the images were recorded. The median of the alumina particles 
is approximately 10 /*m. Due to the recording process, some of the particles look like 
open circles. Although they are all non-overlapping in space, some particles close 
together appear to be connected due to the projective nature of the recording. 

The planar pattern of particles is quite sparse for which reason it might be modeled 
by a thinned marked Poisson process as in Sections [3] and |4j However, a Matern-II- 
type model from Section [4] might be comparatively more promising due to the higher 
attainable intensities, see Remark [T3| 

Our first attempt was to fit a Matern II process for hard spheres as in Remark [10] 
with gamma-distributed radius marks, where for practical reasons the distribution was 
truncated at some high value. The distribution of the weight marks was chosen to be 
the uniform distribution on [0, 1]. Here, three parameters, (A — 315 mm -2 , shape= 6.5, 
rate= 0.00128 mm), were estimated again by the minimum contrast method using the 
pair correlation function. A comparison of the resulting model pair correlation function 
gM and the empirical pair correlation function g of the data shown in Figure [8]indicates 
that this kind of model is not flexible enough already for the second-order behavior 
of the data. This is supported also by the visual inspection of the corresponding 
('Matern') nearest-neighbor distance distribution function and empty space function, 
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Figure 8: Plot of g (dotted), gu (dashed) g t h (solid), left. Function f c for fixed c and 
to, n. 
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Figure 9: Empirical p.d.f. of radii of the data (dotted), Matern (dashed) and n t h 
(solid). 



respectively, in Figure [TO] see also Table [T] for several related deviation tests. 

This motivated us to try modeling with the generalization introduced in Section [4] 
Here, the ansatz is n t h ~ MatII[A, T(a, (3), Uniform[0, l],Po,f c ] with 



f c (r,m,n) 



1 < r < to + n 

exp[— c(r — to — 77.)] otherwise 



(shown in Figure [8| and T(-, •) as gamma distribution. 

Parameter estimation by the minimum contrast method (11) with the pair correla- 
tion function yielded estimates A = 335 mm~ 2 , a — 6.3, j3 — 0.00127mm, po = 0.96 
and c = 119 mm -1 . The resulting characteristics like the pair correlation function <?th 
shown in Figures [8| to [To] indicate a much better fit than the model of the first attempt. 
This is also supported by the corresponding deviation tests, see Table [l] each based 
on fc = 99 simulations of the fitted model. 



6 Conclusions and outlook 

In this paper, we have examined a new class of point processes generalizing the Matern I 
and II hard-core point processes as well as the independent thinning approach. Clearly, 
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Figure 10: Left, empirical c.d.f. of the nearest-neighbor distance of the data (dotted), 
Matern (dashed) and II th (solid). Right, same for empirical empty space 
function. 





pcf 


p.d.f radii 


nearest-n. 


empty sp. 


Matern 


0.00 


0.13 


0.01 


0.00 


n th 


0.16 


0.17 


0.07 


0.06 



Table 1: p- values for the deviation tests of the models 'Matern' and n t h using the 
pair correlation function (pcf) , the probability density function of the radius 
marks after thinning (p.d.f radii), the nearest neighbor distance distribution 
function (nearest-n.) and the empty space function (empty sp.), each based 
on 99 simulations. 
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the proposed new model is not suited for very dense and structured packings of parti- 
cles, since, like for the Mat em processes, there is a relatively small upper bound for the 
intensity for any given / not vanishing almost everywhere. However, it provides a flex- 
ible and simple to fit model in the class of dependent point process models. Unlike for 
the also very flexible Gibbs point processes the simple mathematical structure allows 
for a simple and straight forward simulation and a direct computation of structure 
functions of the resulting point process, also simplifying the application of standard 
fitting procedures. The approach allows for a interpretable descriptions of interac- 
tions, like interaction of shaped objects or non-deterministic death from competition. 
In our further research we also found that this class corresponds to distributions de- 
rived for snapshot point patterns in moving particle systems models we developed for 
the alumina particles. However, this relation has to be discussed in a separate article 
introducing these moving particle systems modeling. 

Thus we think that this new class is worth considering for modeling various real 
world point patterns and systems of particle centers. 

7 Proofs 



7.1 Proof of Theorem M 

The intensity measure 



Waagepetersen 2003 Stoyan et al.[ 1995 a t h of ^th satisfies 



Daley and Vere- Jones 1988 Illian et al. 2008 M0ller and 



a t h(d(x, m)) = po q m a*(d(x, to)) 



(12) 



where a^(d(x, m)) = Xdxfi(dm) is the intensity measure of "J and q rn is the proba- 
bility that a primary point with mark to in ^ is retained as secondary event in 
when thinning would be restricted to pairwise interaction, i.e., (JlJ would be applied 
with po = 1 . 

Using Palm theory [Daley and Vere- Jones 1988 Illian et al.| |2008 Stoyan et al. 
1995] and stationarity of VP, q m is the probability that under the reduced Palm distri 



bution P, n of 'J, where o denotes the origin, the point (o, to) is not deleted by any 
other point. Since \& is a marked Poisson point process this is, due to the Slivnyak- 
Mecke theorem Illian et al.l 20081 M0ller and Waagepetersen 2003 , equivalent to the 



probability that the point (o, to) is not deleted by any point from when the same 
thinning rule is applied. 

Let Q m be the marked point process which consists of all points from W causing a 
deletion of (o, to). Then q m is simply the probability that Q m has no points. Obviously, 
Q m is obtained by independent thinning of "J, i.e. 



^2 P(x,m,l)S {xJ) , 



where /3(x, to, I) is Bernoulli-distributed with parameter /(||x||, to, I) and fyz,/) denotes 
the Dirac measure centered on (x, I). Hence, Q m is an inhomogeneous marked Poisson 



1G 



process Illian et aT| 2008 Section 6.1] with intensity measure 

a Qm 0) = /(II 21 II) m ; a*(d(x, i) = /(||a;||, m, I) Xdx /i(dm). 
Since q rn is the void probability of the Poisson process Q m this implies 
q m = P(Qrn(R rf x R) = 0) = cxp (-a Qm CR d x R)) 
= exp I -A / / 

exp ( —Xdbd 

R ^0 



/(||x|| , m, Z) da; /i(dZ) 

/(r,m,0r d_1 dr/i(d0 



(13) 



using polar coordinates in the last step. Hence, due to Equation ( 12 ), the point process 
of unmarked points of ^th has intensity 



Ath = Ap / q m fi{dm) 
= Xp J exp (-\db d 



/(r,m,0r d - 1 dr A t(d0 n(dm). 



R JO 



7.2 Proof of Theorem M 

Let K m ,n{r) be the probability that two points in 'J with mark to and n a distance r 



factorial moment measure Daley and Vere- Jones] 


1988 


Illian et al. 


Waagepetersen 


2003 


Stoyan et al. 


1995 


a th °* ^ th satisfies 



a lh(dO, m , y,n))= Po 2 K m , n (\\x - y\\) o$(d(x, to, y, n)), (14) 

(2) 

where the second-order factorial moment measure of W factorizes to 

(d(x, to, y, n)) — a*(d(x, to)) a*(d(y, n)) — A 2 dx [i(dm) dy n(dn) (15) 

since * is a Poisson process [see Daley and Vere-Jones[ 1988] . Using again Palm theory 
and the Slivnyak-Mecke theorem, K m!n (r) equals the probability that the two points 
(o, to) and (z, n), \\z\\ — r, do not delete each other and are non of them is deleted by 
any point from ^ according to the thinning rule Q with po — 1. Since both events 
are independent and the probability of the first event is (1 — f(r, to, n)) 2 , it follows 



S»W = (1 - f(r, m, n)) 2 P(W^, m ,„(R d x R) = 0), 



(16) 



where W z ^ m ^ n is the marked point process which consists of all points from "J causing 
a deletion of (o,to) or [z, n). Due to independent thinning, i.e., 

W z 

,m,n 

(a,/)e* 



17 



where j r (x,a,b) is Bernoulli-distributed with parameter f(\\x — r ■ v\\,a, b), \\v\\ = 1, 
Wz,m,ra is an inhomogeneous marked Poisson process with intensity measure 

«Wi, TO , n (d(x,i)) = [/(||ar||,m,i) + /(||x - r ■ v\\,n,l) - f(\\x\\,m,l)f(\\x - r ■ v\\,n,l)] 
x Ada; n(dl). 



According to Equation ([16]), this yields 

SnH = (1 - f(r,m,n)) 2 exp (-a Wzmn (]R d x R)) 

= (1- f(r,m,n)) 2 exp(-X [ f f(\\x\\ , m, I) dx /i(d/) 

x exp ^— A 

x exp ^A 

= (1 - f(r,m,n))' 2 q m q n exp I A y /(•.///./)• /(■.//./ )l /•)//(<!/ ) 

using Equation ( 13 1 and the radial convolution Equation ([3]) in the last step. Ab- 
breviating the last factor of the product in Equation (|17| by q m ,n(f) and combining 



f(\\x\\,n,l)dxiJ,(dl) 

Z)/(||ac - r ■ v\\,n,l) dx(i(dl) 



'R JR d 

^2, 



(17) 



equations Equation (14 1 and Equation (15 1, the second-order product density of the 
point process ^th of unmarked points of is 



0th(r) = X 2 Po' 



(l-/(r, m, n)) 2 q m q n g m ,„(r)^(dm)f((dn). 



R JK 



Due to gth(r) — Qth(r) / A^ h for r > [illian et aT 2008 and Theorem lathis yields the 
asserted form of the pair correlation function <7th> i. e., in particular, po cancels out. 



7.3 Proof of Theorem 



Basically the idea of the proof is the same as in Section 7.1 Here, let q m {w) be 
the probability that the point (o, m, w) is not deleted by any point from II when the 
thinning rule Q with pq = 1 is applied. Then po q m (w) is the density of the intensity 
measure of n t h with respect to the intensity measure an of II, aji(d(x, m, w)) — 
Xdx v m {dw) fj,(dm), and q m (w) equals the probability that the marked point process 
Qm,w consisting of all points from II causing a deletion of (o, m, w) is empty. Since 
Qm,w is a Poisson process with intensity measure l{v < u>}/(||ie||, m, /) au(d(x,l,v)) 
this yields 

q m (w) = P(Q m , w (R d x R x R) = 0) 

l{v < w}f(\\x\\,m, I) dxvi(dv) n(dl) 

F lyi (w)/(r,m,/)r^ 1 dr/i(d/) ) (18) 



exp —A 



'R JK JR d 

exp ( -Xdb d / / 
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and, finally, 



\th = Ap / / q m (w) v m (dw) fi(dm). 

JK JR. 



7.4 Proof of Theorem [14] 

The main arguments of the proof of Theorem [9] in Section |7.2| can be carried over. 
Let po 2 K m , n (w,t 1 r) be the probability that two points in II with marks m and n 
and weights w and t a distance r apart are both retained in n t h- Then K m ^ n (w, t,r) 
equals the probability that (a) the two points (p,m,w) and (z,n,t), \\z\\ = r, do not 
delete each other and (b) non of them is deleted by any point from II according to the 
thinning rule ^ with p — 1. Again, both events are independent, and the probability 
of event (a) is 

1 - [TP(A) + F(B) — P(A n B)] = 1 — f(r, m, n) + l{t = w}[f(r, m, n) 2 - f(r, m, ri)] 

since the probabilities that (A) (z,n, t) deletes (o,m,w), that (B) (o, m, w) deletes 
(2, n, i), and that both delete each other are l{t > w}f(r,m,n), l{t < w}f(r,m,n), 
and l{t = w}f(r,m,n) 2 , respectively. The probability of event (b) is the probability 
that the Poisson process W-^m.n.tu,* consisting of all points of II causing a deletion of 
(o, m, w) or (z,n, i) is empty. Hence, using q m (w) from (18) as shorthand, it equals 

P(W 2 ;, m , n , w ,t(M d x R x R) = 0) = q m (w)q n (t) q m , n (w,t,r), 

where 



q m ,n{w,t,r) = exp A / / F Vl (min{w,t}) f(\\x\\,m,l)f(\\x — r ■ v\\,n,l) dx^(dZ) 

since W T z ,m,n.w,t has intensity measure 

[l{u < w}f(\\x\\,m,l) + l{u < t}/(||a;|| > n,0 

- l{u < w}f(\\x\\,m,l)l{u < t}f(\\x\\,n,l)] a n (d(x , I , u)) 

and 

/ l{u < w}l{u < t}vi(du) = F y; (min{w,t}). 
Jb. 

Therefore, the second-order product density of is 

Qth(r) = \ 2 po 2 / / / K m ,n(w,t,r)v m (dw)v n (dt) /J,(dm) (J,(dri) 
Jb. Jr. Jr. Jr. 

— X 2 po 2 / / (1 — fir, m, n)) I r (m, n) /i(dm) fi(dn) 
Jr Jr 

where 

I r (m,n) = / / q m (w) q n (t) q rn , n {w,t,r) v m {dw) v n {dt) . 
Jr Jr 

Note that the summand l{t = w}(f(r,m,n) 2 — /(r, m, n)) in the first factor of the 
integrand has been left out since its integral vanishes due to the assumed continuity 
of the distributions v m and v„ . 
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